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Introduction 


Mathematical modelling of reality usually results in an equation or a set of 
equations. These equations can have a variety of forms, and the nature of the 
solutions to the equations depends on the nature of the equations. 


You should already be familiar with ordinary algebraic equations, whose solutions 
are numbers. For instance, the algebraic equation 
x?-2=0 


has a solution 


x =,/2=1.414.... 


In Unit 1, you studied equations of a different type called recurrence relations. 
These are equations whose solutions are sequences. For example, we saw that the 
sequence whose nth term u, is defined by 


u, = 2" 
is a solution of the recurrence relation 
u,, = 2u, (se 1, 2, 353.3) 


The equations we study in this unit are of a different kind again; they are called 
differential equations. Such equations are widely used in mathematical models. You 
will see examples in many units of this course, starting with Units 3 and 4. An 
example of a differential equation is 

dy 


gee 1 
Fe a (1) 


You will see that it contains two variables x and y. The problem posed by a 
differential equation is not to find numerical values such as Je for the variables in 
it but to find a function relating these variables. One way of describing such a 
function is to give a formula expressing y in terms of x, for example the formula 


y=e-x-1l (2) 


describes one such function. In fact the particular function described by Formula 
(2) happens to be a solution of Equation (1); for if x and y are related by Formula 
(2) then the left-hand side of Equation (1) becomes 
dy _ 
dx 


and the right-hand side becomes 


d 
ax ° —x-l)=e-1 


X+ty=xt+(F—-—x-1l=e-1. 


That is, Formula (2) makes the two sides of Equation (1) equal and so the formula 
gives a solution of the differential equation. 


By the end of this unit you will know how this solution was obtained, and you 
will be able to find all the other solutions of the differential equation as well. 


There is an alternative notation for differential equations and their solutions. —— 
I ; ; This is the notation used in 
nstead of using two variables, x and y say, this alternative notation uses just one M101 and MS283. 

variable, x say, together with a symbol for the function the differential equation 

informs us about. For example, if this function were denoted by U, then Equation 

(1) would take the form 


dU (x) ere 
dx 


or 
U'(x) =x + U(x) 
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and Formula (2) would be written 
U(x) =e*—-x-1 
or 
U:xrroe—x-1. 
In this course, however, following the usual practice of applied mathematicians, 


scientists and technologists we shall not normally introduce a symbol for the 
function. 


To translate from one notation to the other, use the following ‘dictionary’ 


this course function notation 
y U(x) 

dy dU(x) 

7 aE or U'(x) 


Study guide 


There are two main methods for solving differential equations: analytical methods 
which provide the solution function as a formula like Formula (2) above, and 
numerical methods which provide the solution function in the form of a table. In 
addition, there are graphical methods which provide information about the graph 
of the solution function. 


This unit consists of two parts which can be studied independently. One part, 
consisting of Sections 1, 2 and 6, is concerned with numerical and graphical 
methods. It begins with a television programme on some basic ideas about 
differential equations, including a discussion of some graphical techniques; it then 
proceeds, in Section 2, to describe a simple numerical method. Section 6 describes 
the work in your first computer terminal visit, part of which is designed to help 
you learn to use the numerical method described in Section 2. Be sure to plan this 
terminal visit in good time: there may be a deadline after which you cannot use 
the terminal. For up-to-date information about this consult the Course Guide and 
stop press notices. 


The other part of the unit, consisting of Sections 3 and 4, describes some 
analytical methods for solving differential equations. You can do this part of the 
unit either before or after the television programme and computer terminal visit. 
There is a tape commentary which forms part of Section 4. 


Section 5 contains exercises to enable you to revise and consolidate what you have 
learned from the rest of the unit. There is no essentially new material in this 
section. 


The order in which you study the sections of this unit depends on how the 
television broadcast schedule fits in with your study schedule and also on when 
you can arrange your visit to the computer terminal. If you are ready to start 
reading the unit before the television programme is due to be broadcast, begin by 
studying Sections 3 and 4, returning to Sections 1 and 2 after watching the 
programme. Otherwise start by watching the television programme and then study 
the unit in the order it is written. Whichever part of the unit you read first make 
sure you read through the preliminary work at the beginning of Section 1 before 
watching the television programme. Before going on your computer terminal visit 
you must have studied Section 6. This can be done any time after studying 
Sections 1 and 2. 
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1 Direction fields (Television Section) 


1.1 Preliminary work 


This section describes an approach to differential equations starting from the 
concept of direction field. This concept is explained in the television programme, 
but the following exercise will help you to understand the basic ideas more easily. 


Exercise 1 

On the graph paper provided in Figure 1 draw short line segments 

(i) through (0,0) with slope 1 

(ii) through (3,4) with slope 0 

(iii) through (, —4) with slope —# 

[Solution on p. 43 | 

The television programme looks at the practical problem of finding curved mirrors 
which focus a parallel beam of light towards a particular point (see Figure 4). 
Using this problem as an example the programme shows how any direction field 
leads, in a natural way, to a set of curves (called trajectories). The method used to 
find these curves is illustrated by various examples. Exercises 2 and 3 ask you to 
check some of the results used in these examples. 


Exercise 2 
(i) Given that y = 4x? + 5, check that 
dy 
ae x 
(ii) Given that y = 4x* + C, where C is any constant, check that 
dy 
ao x 


(iii) Given that y = Ce* — x — 1, where C is any constant, check that 


dy 
a 


[Solution on p. 43 ] 


Exercise 3 
(i) Given that y = (R* — x”) '/*, where R is any constant, check that 


‘ye 

a 
To keep the calculation under control I suggest the following systematic 
procedure. 


d 
(a) Calculate = from the given formula for y. 
x 


(b) Calculate = from the given formula for y. 


(c) Check that the results of (a) and (b) are the same. 
(ii) Given that y = = a. with C > 0, use a similar systematic procedure to check that 


oy yt ey 


dx x 
[Solution on p. 43] 


During the programme we shall need to find the orientation of a reflecting surface 
which will deflect a ray of light towards a particular point. To get a formula for 
this orientation we need the law of reflection (see Figure 2). This is an empirical 
law which tells us that the normal (i.e. the perpendicular) to the reflector bisects 
the angle between the incident and reflected rays at the point of reflection. That is 
to say, the angles a and £ in Figure 2 are equal. Since « + 6 and f + ¢ are also 
equal (both being right angles) it follows that 


O= 9g. 


This relation is used as part of a geometrical derivation in the programme. 


reflected ray 


y = (R? — x”)"’? is the 
equation of a semi-circle of 
radius R. 


incident ray 


normal .. 


point of reflect 


reflector 


Figure 2 
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In order to follow this derivation you also need to know the following two 
geometrical theorems (for how they are used, see Figure 9). 


O QO 9 [sp 


Figure 3(b). If the angles w and 6 are 
equal, then the lengths of sides OQ and 
OP are equal. A triangle such as OOP 
with two equal sides is said to be 
isosceles. 


Figure 3(a). If the lines OQ and PR are 
parallel, then the angles W and 6@ are 
equal. 


Now view the television programme ‘Direction fields and families of curves’. 


1.2 ‘Direction fields and families of curves’ 


In this programme we look at the problem of focusing a parallel beam of light to TV? 
one point called the focus: we want to formulate a mathematical equation that 
describes a curved reflector which will do this (see Figure 4). 


parallel beam 


ae reflector 
Nel” 


To analyse the problem, we look at one ray of light selected from the beam and 
show that if a small reflecting surface is placed at a point, then there is only one 


orientation of this surface for which the ray is deflected towards the focus (see 
Figure 5). 


Figure 4 


Figure 5 


By considering other points in the plane, we obtain a whole set of directions. 

Every point in the plane (except for the focus) has an associated direction which is You can also think of a_ 
a 5 direction field as a function 

unique and we can determine it by experiment. Such an association of each point 


Se ee which maps points in the 
with a unique direction is called a direction field (see Figure 6). plane to slopes. 


ee et ee. SS. 
ee eee 
‘oe eS fof 2s st 
—. = — 7 
NOUN ew OF 
i a a 
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At the beginning of the programme, we showed three curves that had the required 
property of focusing a parallel beam of light onto one point. By superimposing 
one of these curves on the direction field (see Figure 7), we immediately notice 
that the slope of the curve at any point coincides with the slope given by the 
direction field. To illustrate how this observation can be used to find equations for 
the curves we first consider a simpler problem. 


We consider the problem of finding equations for the curves that fit the direction 
field given by the formula 


slope at (x,y) = ee 
b 2 


A diagram of this field is shown in Figure 8. 


te ee pe eS 
ed i 
ae er wis SK NNN 
oy See fee 
a. NS el we Sf 
a el oe Oe 
ee ee ee oe 


Figure 8 


Question 
Sketch on Figure 8 a few curves that fit the direction field. 


This is done by moving your pencil so that it is always travelling in the local 
direction of the field. 


The sketch suggests that the curves fitting the direction field are circles centred at 
the origin, and the computer graphics confirm this; but we would like to be able 
to check this visual judgement by exact calculation. To check that the circles do 
fit, we calculate the slopes of the tangents to the circles. The equation of a circle of 
radius R, centred at the origin,is x7 + y? = R* and so expressing y as a function of 
x we obtain (for the part of the circle where y > 0) 


y = (R? — x?yt2, 
Differentiating, we obtain for the slope of the tangent at (x, y) 
a —2x x 
aa Se - ar] — ae see Exercise 3(1). 
So we have checked that the family of circles do fit the direction field whose slope 


eee 
at (x,y) is —-—. 
y 


Now we return to our original problem of focusing a parallel beam of light. The 
example above suggests that it would be helpful to have a formula for the 
direction field. The slope of the direction field at (x,y) is the slope of the mirror 
which reflects the vertical ray of light through the focus O (see Figure 9). By 
considering the triangle PRQ we see that 


We are using the convention 
PR ee PS + SR ae = at that, for example, PR is the 


OR ~. OR x length of the line from P to R. 


slope at (x,y) = 
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Since OORS is a rectangle this becomes 


y+0OQ 
a 


slope at (x,y) = (1) 


Now consider the angles. Since OQ and PS are parallel we have (see the theorem 
at the end of Subsection 1.1) 


yy = ¢. 

Also, by the law of reflection we have 
Ge. 

It follows that 
y= @. 


Thus two angles in triangle OPQ are equal, and so the corresponding sides are 
equal: 


OO = OP 
= ./x* + y’ (by Pythagoras). 


Substituting this into Formula (1) we see that the formula for the direction field is 


y+ ./x?+'y* 


xX 


slope at (x,y) = 


What condition must be satisfied by a curve that fits this direction field? Suppose 
the equation of the curve is 


y=f (x) 


where f is a function we do not know yet. The slope of a tangent to the curve is 
dy/dx, the derivative of f(x). If the curve is to fit the direction field this slope must 
match the slope of the direction field, that is, we must have 


ope Ey (2) 


dx x 


An equation like this, connecting the derivative = of some unknown function to 
x 


the variables x and y which are related by it, is called a (first order) differential 
equation. 


There are many different curves that fit the direction field, and each one of them 
corresponds to a different function whose derivative satisfies the differential 
equation. These functions are called solutions of the differential equation. 


Our original reflector problem has now been reduced to the problem of finding 
solutions of Differential Equation (2). At first sight, finding the solutions of a 
differential equation may seem like a formidable task. Fortunately it is often 
possible to find a method which gives all the solutions in terms of a single 
formula. To see how this can be done, we again digress from the reflector problem 
and consider a simpler differential equation 


| 
— 


The solutions of this differential equation can be found by integrating both sides 
with respect to x to give the formula, 


y=3x°+C 


where C is an arbitrary constant of integration. Each numerical value for C gives a 
different function relating y to x, and as you checked in Exercise 2 each of these 
functions satisfies (i.e. is a solution of) the differential equation. Each of these 
solutions corresponds to a different curve fitting the direction field whose slope at 
(x,y) is x (see Figure 10). 


Figure 9 

y 
1e ks oer 
\X+ 4 
ee oe ale 
WA + 47 


SEY 


yok *® ve 
Figure 10 
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A more complicated example is given by the differential equation 


dy 
S eae 


The solutions of this differential equation can be found by the method described in 
Section 4 of this unit. They are described by the formula 


y=Ce—-x—-l. 


Since the process used to obtain this formula is more complicated than a direct 
integration, the arbitrary constant C is this time not just added on to the 
expression for y but is embedded in it in a more complicated way. 


Returning again to the reflector problem, we want to solve the differential 
equation 


ly ye ay 
dx x 


This is not one of the standard types considered in this unit, but it happens that 
the equation can be solved by introducing a new variable u, defined by 


k= RY. (3) 
(The details of the following calculation are not in the television programme and 
are not essential to the course). Squaring both sides of Equation (3) and then 
differentiating with respect to x gives 
du dy 
2u— =2 2y— 
‘ie ae 


x? 2 
= 2x +2yP4VE HY 
x? + y? + y,./x? + y? 


> 


Xx 


(from the differential equation) 


=2 


2 
=o a (by definition of u) 


so that 
du is, 
ee 
fx ty ty 
x 
dy 
ae (from the differential equation) 
x 


It follows by integrating both sides with respect to x that 
u=y+C 


where C is a constant of integration, which must be positive since 
u = ./x* + y* > y. Squaring both sides and then using the definition of u gives 


x? +y? =y? + 2Cy+ C? 
which simplifies to 


yo 20 y ee C3 
a & 
1.€. a 


This is the formula for the solution of the reflector problem used in the 
programme. 
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Summary of Section 1 
The main points made in the programme are these: 
1. A direction field specifies a slope at each point in the (x,y) plane. 


2. A first order differential equation is an equation connecting the derivative 
dy/dx of some unknown function to the variables x and y which are related by 
that function. 


If y = f(x) is the equation of a curve that smoothly fits the direction field, then 
the function f must satisfy the differential equation obtained by equating dy/dx 
to the expression for the slope of the direction field. Any function (that is, any 
relation giving y in terms of x) which satisfies the differential equation is called 
a soluticn of the differential equation. There is an infinite family of solutions, 
each one corresponding to a different curve that fits the direction field. 


3. If an exact formula can be found for the solutions of the differential equation 
then this formula will contain an arbitrary constant of integration, called C in 
this programme. Different numerical values for C give different solutions 
and therefore label different curves fitting the direction field. 


4. Whether or not it is possible to find a formula for the solution of a differential 
equation, the direction field is a useful source of information about its 
solution; in particular it can be used to sketch solution curves and so get an 
idea of their qualitative behaviour. 


Exercise 4 
Draw the direction field whose slope is 2 at every point in the plane. Can you deduce the 


family of curves associated with this field? Check your result by differentiation. 
(Solution on p. 43 | 


Exercise 5 


On p. 8 we checked that the upper half of the circles matched the direction field. Can you 


do the same for the lower half? 
(Solution on p. 43 | 


Exercise 6 


What would be the differential equations for the trajectories of the direction fields defined 

by 

(i) Slope at (x,y) = —x/y 

(ii) Exercise 4 

(iii) The direction at any point P (other than the origin) is the same as that of the line OP 
where O is the origin. 

(Solution on p.43 ] 


Exercise 7 


For the direction field shown in Figure 11, sketch a few trajectories, including the one that 
passes through the point P. 


Figure 11 ee XN 
a ee io eet 
| le fie Re es 


So eet 
se ow 


‘i 
x 
f 
f 


SN ON 


fee ee 
(Solution on p. 44 | fF og _- 


f 


Exercise 8 


Using the graph paper in Figure 12, draw some slopes belonging to the direction field 
defined by 


slope at (x,y) =1-—y 


These curves are often called 
trajectories. 


These functions are sometimes 
called particular solutions. 


This formula is called the 
general solution of the 
differential equation. 
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including the ones at the points (0,0), (0,5), (0,1), (0, 14), (0,2), (1,0), (1, 1), (1,2). Hence 
sketch enough slopes to give a good idea of how the direction field behaves in the region 
shown. 


Figure 12 


ex aT i + 2s GRP a5 aon flay ae ze 


: j 


H H j 


oe 


! 
‘ 


[Solution on p. 44] 


inti 
2 x 


Exercise 9 


Use your direction field from the preceding exercise to sketch some solution curves of the 
differential equation : 


including the ones passing through the points (0,0), (0,1) and (0, 2). 
[Solution on p. 44] 


2 A numerical method 


2.1 A geometrical description of the method 


So far the only systematic method we have for solving differential equations is to 
draw the direction field and sketch its trajectories. In the rest of this unit we 
develop some other methods for solving differential equations. 


For certain types of differential equation there are exact methods of doing this 
which give the solution as a formula such as y = e* — x — 1. Some of these 
methods are considered in Sections 3 and 4; but there are many differential 
equations for which these exact methods do not work and an approximate 
(numerical) method of solution must be used. Most numerical methods yield the 
solution not as a formula but as a table of values, similar to a table of logarithms. 


In this section we shall study the simplest numerical method for differential 
equations. The method is known as Euler’s method, after Leonhard Euler 
(pronounced Oiler), who lived from 1707 to 1783 and was one of the most prolific 
mathematicians of all time (his complete works fill some 60 to 80 volumes). Euler 
devised this method for computing the orbit of the moon; it was used to compute 
lunar tables for the British Admiralty, Euler being voted £300 for his method, 
while the person who did the calculations received £5000. Euler’s method is not 
used much now but the more efficient methods that have been developed since are 
based on similar principles. Details of some of these more efficient methods are 
given in the unit on numerical solutions of differential equations. 


—— 


— — = Be 
ee ag Sig 
- ~~ ~~ 
ae sa 
> 
Z a _ 
we 
3 ~~ 
— ee 
A =e te 
oe ~~ 
oe 
, a 


MST204 2.2 13 


The idea behind Euler’s method is to approximate the trajectories of the direction 
field associated with the differential equation by segments of straight lines joined 
end to end to form an open polygon as shown in Figure 1. 


Before discussing how to construct this polygonal approximation we must first 
specify which of the infinite family of trajectories we wish to approximate. Let us 
suppose that the trajectory has been specified by saying that it passes through a 
certain point P, with co-ordinates (x, y,). Such a specification is called an initial 
condition; it provides a starting point from which the polygonal approximation 
can be built up step by step as follows. The first vertex of the approximating 
polygon is taken to be the point P, itself. The first of the line segments forming 
the polygon starts at P, and is taken to have the same slope as the direction field 
at Po; let us call this slope m, (see Figure 2). This choice makes the first segment a 
tangent at Py to the trajectory we want, and gives it the right slope at P,; but as 
we proceed along this line segment its slope will deviate more and more from the 
slope of the direction field. If this error is not to become too large, we must soon 
make a change of direction. Let us call the point at which we decide to make this 
change of direction P,: it is the next vertex of our polygon (see Figure 3). At this 
point we begin the next line segment of the polygon, giving it the slope of the 
direction field at P,, which we denote by m,. This new segment is continued to 
some point P, (see Figure 4) at which we once again change direction, the next 
segment having slope m, equal to the slope of direction field at P,. Continuing 
this process we obtain the open polygon shown in Figure 1. 


—— 


<_ Soe ~~ = 
a= ~ ~~ 
—_ pecs ~~ ~~ 
= o = ~ 
a 
a 
> S 
a 
~~ 
a4 
“A ~~ 
4 
Figure 3 Figure 4 


To follow the trajectory perfectly we would have to adjust the slope of our 
constructed figure continuously to that of the direction field, but in this 
approximation we adjust it only at the points P,,P,,P,... . 


Exercise 1 


The diagram below shows a direction field whose slope at the point (x, y) is x + y. Using 
the graphical construction just described, draw on the diagram a polygonal approximation 
to the trajectory passing through the origin P,. Choose the points P,, P,,...,P, so that 
their x co-ordinates are 0.2,0.4,...,1.0. There is no need to make any detailed calculations; 
a rough sketch is all that is required. 


oe 


<<. 


Po 
U2 04° 06 - DE. *EO 
— ss, — er 
~~, ~~ —— oo 
Figure 5 


(Solution on p. 44] 
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2.2 A formula for Euler’s method 


If Euler’s method were no more than a graphical construction, it would offer little 
advantage in comparison with sketching solution curves directly from the direction 
field, which is the method we used in Section 1. The point of Euler’s method, 
however, is that it can, with advantage, be carried out numerically rather than 
graphically. Instead of a graph showing the polygonal approximation to the 
trajectory, the numerical approach yields a table listing the x and y co-ordinates of 
the vertices Py, P,,.... By referring to this table we can read off the 

approximation for y corresponding to the value of x we are interested in. 


In order to simplify the calculation of the x co-ordinates we choose the vertices 
Po, P,,... so that their x co-ordinates X9,x,,... are equally spaced. The spacing is 
called the step length and is denoted by h. For example, in Exercise 1 the x co- 
ordinates were 0,0.2,0.4,... and so the step length h was 0.2. In general, suppose 
the trajectory we are trying to approximate has been specified by giving the co- 
ordinates (Xo, yo) of the point P,. If the step length is h then the x co-ordinate of 
P, is given by x, = xX, + h. By adding another step length h we obtain the x co- 
ordinate of P,; x, =X, + 2h. Continuing in this way we obtain the general 
formula, 


X, =Xq + rh 
for the x co-ordinate of the vertex P.. 


The y co-ordinates of the vertices are calculated by setting up a recurrence relation 
from which the y co-ordinate y,, of the vertex P,,, can be calculated once the y 
co-ordinate y, of the preceding vertex P, is known. To obtain this recurrence 
relation we use the condition formulated in the preceding subsection that the slope 
of the line segment P,P,, , is to equal m,, the slope of the direction field at P. (see 
Figure 6). 


| 
| 
| 
| 
Ce Ee OE ie 


. eee we KONO 
aes SY XS 


Figure 6 


-E r+] 


From the right-angled triangle shown in Figure 6 the condition for the line 
segment P,.P,,, to have slope m, is 


length of P,,,Q _ = 
lengthofPQ ”™ 


Now the length of P,, ,Q is the difference between the y co-ordinates of the 
vertices P,, , and P,, that is y,, , — y,. The length of P.Q is the difference between 
the x co-ordinates of P,, , and P,; it is equal to the step length h. We can 
therefore write 

Bf r+1 y — 


mM,. 
h r 
Solving for y,,, we obtain the recurrence relation for the y co-ordinates 
Vr+1 ge yr + hm,. 


Since the y co-ordinate of the point P, is known, the y co-ordinates of the 
remaining vertices can be calculated using this recurrence relation. However, 
before we can do this we need a formula for the slope m, of the direction field at 
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the point (x,,y,). This formula is provided by the differential equation we are 
trying to solve: for example if the differential equation is 


dy 
Te 


then the formula for m, is 
m, =X, +y,. 


For a general differential equation of the type considered in this unit, say 


dy 
dx ne m(x, y) 


where m(x, y) can stand for any formula involving one or both of the variables x 


and y, the corresponding formula for m, is 
m, Ba m(x,, y,). 


To sum up we have: 


Procedure 2.2: Euler’s method 


To apply Euler’s method to the differential equation 


dy 
dx 4 m(x, y) 


given the initial condition, y = y, when x = Xo, using step length h; 

1. Use the co-ordinates x,, yo as a starting point. 

2. Calculate the next y co-ordinate using the recurrence relation 
Yui = y, + hm, 

where m, = m(x,, y,). 

3. Calculate the corresponding x co-ordinate using the formula 
X= 3, +28 

(after increasing the value of r by 1). 


4. If more co-ordinates are required return to Step 2. 


The recurrence relation for Euler’s method can also be derived without reference 


to a graphical construction. By the definition of a derivative, the differential 
equation is equivalent to 


lim 2% + h) — y(x) | 
1m —---——_- = 
h>0 h 


m(x, y). 


Hence, for sufficiently small h, we have the approximation 


y(x + h) — y(x) 


; ~ m(x, y). 


15 


This holds for all values of x, in particular for x = x,; so we can replace x by x, in 


the formula and rearrange to obtain 


y(x, + h) = y(x,) + hm(x,, y(x,)). 


With the notation y(x,) = y,, v(x, +h) = y,,, this is the same as the recurrence 


relation in the ‘Euler’s method’ box. 


Example 1 
Let us apply Euler’s method numerically to the problem treated graphically in 


Exercise 1. In that exercise we considered the direction field corresponding to the 


differential equation 
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== x+y 
and constructed a polygonal approximation to the trajectory passing through the 
origin. The step length was 0.2. 


Since our polygonal approximation started at the origin, we take P, to be the 
origin. So the starting point for our calculation is 


Xo = 9, Vo = 9, My = Xo + Yo = 9. 
According to Procedure 2.2, the next y co-ordinate is 
Vy =Yo + 0.2m, = 0 
and the corresponding x co-ordinate is 
X, =X,+1x02=0.2 
which tells us all we need to know about P,: 
x, = 0.2, y, =9, m, =xX,+y, = 0.2. 
Next we obtain the co-ordinates of P,: 
y, =y, + 0.2m, = 0.04 
and ) 
Xo =X,+2x02=04 
so we have 
x, = 0.4, y, = 0.04, m,=X,+ y, = 0.44, 
and so on. 


If you have to do more than a very few steps of such a calculation by hand, it is a 
good idea to lay it out as a table, for example: 


After each value of y,,, has been calculated from the recurrence relation and 
entered in the last column, it is transferred to the y, column in the next row. 


Exercise 2 


Complete the above table and so obtain the coordinates of P,. Check your table by 
comparing it with your graph from Exercise 1. 


[Solution on p. 44 | 


Exercise 3 
What changes would be necessary in the column headings of the above table if the equation 


d 
= = x + y were replaced by 
x 


di 
dx 
For this new equation, use the same step length, h = 0.2, to calculate an approximation to 


y(1) given that y(0) = 1 (the notation y(1) means ‘the value of y when x = 1’ and y(0) means 
‘the value of y when x = 0’). 7 


[Solution on p. 44 | 
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2.3 Choosing the step length: accuracy versus cost 


Whenever we use Euler’s method we have to choose a value for the step length h; 
in the last subsection this choice was made so as to give a convenient illustration 
of the working of the method. When a method such as Euler’s is used for the 
purpose of solving some particular differential equation rather than for illustration, 
however, we have to consider the choice of step length more carefully. Two 
opposing requirements, accuracy and cost, influence this choice, and the decision 
must be a compromise between them. 


The accuracy of Euler’s method usually improves when we take a smaller step, 
size. One way of seeing this is to consider the graphical description of Euler’s 
method used in Subsection 2.1. The smaller the step size, the more often we adjust 
the slope of our polygonal approximation to the trajectory, and so the less the 
approximation deviates from the trajectory (see Figure 7). 


seeeetfipence 


= Ais “ Ssh=0.4 (say) 


_ 
” 
s” 
TA 

Pa 
a” 


Figure 7 


exact solution 


The effect can also be seen by calculation. The next table shows approximate 
solutions, using different values of h, for the problem considered in Exercise 3. 


Given that 


and y(0) = 1, find y(1). 


The table also shows the deviations of these approximations from the exact value 
of y(1), which is e = 2.718282... (you can check that y = e* is a solution of the 
differential equation, that it also satisfies the condition y(0) = 1, and hence that 
y(1) = e* =e). All numbers are given to an accuracy of 6 places of decimals. 


approximation number of 
h to y(1) deviation steps 
0.1 2.593742 0.124539 10 
0.01 2.704814 0.013468 100 
0.001 2.716924 0.001358 1000 


0.0001 2.718146 0.000136 10000 


You will see that the difference between the approximation to y(1) and the exact 
value decreases roughly in proportion to h. Indeed, it can be shown that by 
making h small enough we can make the error as small as we please, in other 
words, the error approaches the limit zero as h approaches zero. This is true of the 
Euler approximation for any reasonable differential equation, not just the one 
considered in the example. 


A word of caution is necessary here. These remarks about the error are only valid 
if we do the arithmetic to enough decimal places; if we are working with a 
calculator or a computer then the number of decimal places we can use is 
restricted and rounding errors will be introduced. After a certain point any 
increase in accuracy brought about by reducing the size of h will be swamped by 
these rounding errors. 
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Rounding errors are not the only problem. Before concluding that h should always 
be taken very small we must also consider the cost of this additional accuracy. The 
last column of the table illustrates how the number of steps goes up in inverse 
proportion to the step size (the general formula for number of steps is (b — a)/h 
where a and 5 are the initial and final values of x). Since the error in Euler’s 
method is approximately proportional to step size it follows that for this method a 
10-fold improvement in accuracy is paid for by a 10-fold increase in the number of 
steps required. For the example above, an accuracy of one part in a million could 
require something like a million steps (and each step would have to be calculated 
to a very high accuracy—perhaps one part in 10'). Even with a high speed 
computer this would not be a sensible way of going about things; instead one 
would use one of the more efficient methods described in the unit on numerical 
solutions of differential equations. Thus Euler’s method is not suitable for high 
accuracy work. Its virtue is rather in its simplicity and its clear illustration of the 
basic principles of the numerical solution of differential equations. 


Exercise 4 


Consider the problem of calculating y(2) (ie. the value of y when x = 2) given that wii} 2, 
for the differential equation 


dy = 
— 


using Euler’s method, for various step sizes. (This equation can be solved exactly to give 
y = 4x’ + 3, but the purpose of the exercise is to test Euler’s method on a very simple 
example.) 


(i) If the calculation is to contain n steps, what is the step size? 


(i1) Calculate the approximation to y(2) for the step sizes that give (a) 1 step, (b) 2 steps, 
(c) 4 steps. 


(iii) Calculate y(2) exactly and hence find the errors in cases (a), (b), (c) above. 


(iv) What do the results of (iii) suggest about the dependence of error on h for this 
problem? 


(v) Assuming your hypothesis in (iv) is correct, how many steps would be required to 
obtain an accuracy of 10° * using Euler’s method? 


[Solution on p. 44] 


Summary of Section 2 


Given a differential equation 


d 
= m(x,y) 


dx 


and the condition that y(xo) = yo, where x, and y, are specified, then Euler’s 
method can be used to find y(b), where b is a given value of x, as follows. 


1. Choose h = (b — x,)/n where n is the number of steps to be used. 
2. Forr=0,1,...,n — 1 use the recurrence relation 
Vr+1 = yp + hm(x,., y,) 
where x, = x) + rh. Then x, = b and y, ~ y(b). 


The error in the approximation to y(b) is roughly proportional to h, and therefore 
inversely proportional to the cost (ie. number of steps) of the calculation. 


3 Direct integration and separation of 
variables 


3.0 Introduction 


This section and the one that follows can be studied either before or after Sections 
1 and 2. They deal with methods of finding exact solutions, that is solutions 
expressible as a formula, for certain types of differential equation. 
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A feature common to all these exact methods is that at some stage they all involve 
carrying out an integration. We shall see that some types of equation have to be 
rearranged before this integration can be carried out but we begin by considering 
some equations which can be solved by carrying out the integration directly. 


3.1 Direct integration I: general and particular solutions 


Here is one example of a differential equation which can be solved exactly: 


dy 
ao = pia 


Before proceeding, see if you can find a solution. Write it in the box below. 


As we have already noted in the introduction to this unit, a solution of a 
differential equation is a function; so we want the box to contain a function of x, 
that is, an expression depending only on x. The condition this expression must 
satisfy is given by the differential equation: it is that differentiation of the 
expression in the box (with respect to x) must give the expression 2x. In other 
words, the box must contain a primitive of the function 2x (a primitive of a 
function f is any function F whose derivative is f). As you probably know one 
primitive of the function 2x is the function x”, and so you could have written x? in 
the box, giving y = x? as a solution of the differential equation. 


This is not the only solution, however; you may recall from an earlier course that 
adding any constant to a primitive gives another primitive. So, for example, the 
function defined by 


y=x*45 


is another solution of the differential equation. (You should check this last 
statement.) In fact, any function defined by 


y= x” +, 


where C is a constant, is a solution of the differential equation. The two solutions 
just considered are the cases C = 0 and C = 5 respectively, but any other 
numerical value for C also gives a solution. 


The method of solution we are using here can be applied to any differential 
equation of the form 


dy 
Pe f (x) 

where f is some given function (for example the function 2x in the above example). 
The important thing is that the right-hand side must not involve y. All you have 
to do in such cases is to integrate the function f. The solutions of the differential 
equation are then given by 


vo [reeax + ‘S 


where { J (x) dx is one of the primitives of the function f, and C is an arbitrary 


constant. 


The appearance of the constant in this solution is a feature common to all 
differential equations; it expresses the fact that there is an infinite collection of 
functions which satisfy the differential equation, one for each value of the constant 
C. The collection of all the solutions of a differential equation is called the family 
of solutions. 


Just as in the case of recurrence relations, we distinguish between particular and 
general solutions of a differential equation. For the differential equation 

dy/dx = 2x studied above, the function defined by y = x? is a particular solution, 
and the function defined by y = x? + 5 is another. But a formula such as 
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y = x* + C, containing an arbitrary constant, is called the general solution. The 
general solution describes not just one solution but the whole family, from which 
we can pick out any particular solution by giving the arbitrary constant a 
particular value. 


Depending on the problem giving rise to the differential equation, we may require 
either the general or a particular solution. If a particular solution is required, then 
the problem will contain some additional piece of information by means of which 
this solution can be picked out from the family of all solutions. For example, 
suppose our problem 1s: 


Find a solution of the equation 
dy 
oe 
such that y(3) = 4. 
The condition y(3) = 4, which means ‘when x = 3, y = 4, is the additional piece of 
information we need in order to pick out the desired particular solution. The way 


to use this information is first to obtain the general solution of the differential 
equation, which we already know to be 


ax 


yox +¢, 


and then use the extra condition y(3) = 4 to give an algebraic equation for C. 
Putting x = 3 in the general solution we obtain 


y=94+C 

and so the condition y(3) = 4 gives 
4=9+4C. 

Solving for C gives C = —5 and so the required particular solution is 
y=x? —§. 

Exercise 1 


(i) Find the general solution of the differential equation 


and the particular solution for which y(1) = S. 
(ii) Find the general solution of the differential equation 
d 
ial a e2* 
_ dx 
and the particular solution for which y(0) = 2. 


(iii) Find the general solution of the differential equation 
dy 
— = asin bx 
dx 


(where a and b are given constants, with b # 0) and the particular solution for which 
y(0) = 0. 


[Solutions on p. 45 | 


3.2 Direct integration II 


All the differential equations we have considered so far in this section have the 
special form 


dy 

=. aS (x) 

_ where f is some given function. Although the derivative of y appears in the 
equation (that is why it is called a differential equation) y itself does not appear. In 
consequence we were able to solve the equation by integrating the function f. We 
can think of this in a slightly different way: in effect, we integrated both sides of 
the equation directly with respect to the same variable x, obtaining 
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( 
\ Sore \J (x)dx +C 3 N.B. Always integrate both 
sides with respect to the same 
variable. 


where y is a primitive of = just as { f (x) dx is a primitive of f; (You might think 


that when integrating both sides of an equation we should also add a constant to 
the left-hand side. However, adding two constants, one to each side of an equation, 
is equivalent to adding their difference, C say, to the right-hand side.) 


Let us now see whether any more general differential equations, containing y as 


well as its derivative, can be solved by integrating both sides of the equation. This equation is equivalent to 
Consider for example, the equation the equation 
dy dy Xx 
2y— = —2x. * Sateen ee 
Y ax dx y 


ae discussed in the television 
We know that the right-hand side can be expressed as the derivative of — x? with programme. 


respect to x. The remarkable thing about this equation is that the left-hand side, 
even though it contains y which is related to x in an unknown way, can be 
expressed as the derivative (with respect to x) of a simple expression. Before you 
read on, try to find a simple expression whose derivative (with respect to x) is 
equal to the left-hand side. (Hint: use the chain rule.) 


The required expression is y’; for if we differentiate y*, using the chain rule, we get 
dy dy 

—— = 2y—. 

dx dx 

This formula is true whatever functional relation connects y and x, so it does not 
matter that we do not know what this functional relation is. 


@5) yee willy”) 
mae eee 


The above remarks show that the differential equation can be written in the form 


d(y*) _ d(—x’) 
eee eae 


It follows that y? and — x? are both primitives of the same function and must 
therefore be equal or differ by a constant. We can therefore integrate both sides of 
the equation to obtain 


y= =x +C 


where C is an arbitrary constant. This last formula contains the general solution of 
our differential equation. 


In practice, it is helpful to lay out calculations of this kind as follows. Starting 
from the original differential equation 


dy 
Y dx < 
we integrate both sides with respect to x, obtaining two primitives which must 


either be equal or differ by a constant. Using indefinite integral notation for these 
primitives, the formula expressing this fact is 


Z —2x 


Remember: there is no need 


fpr) ax= [(-290a + C. to include a constant on the 
x 


left-hand side, one constant is 


Even though we do not know yet how x and y are related, and therefore do not Soe 


know what expressions in x the symbols y and dy/dx stand for, we can use the 
formula for integration by substitution to simplify the left-hand side. In this way 
the equation becomes 


[ava = | -2xax + C. 
Writing this formula in the 


Now straightforward integration gives form y* + x? = C shows that 
> the family of solutions can be 
je ee ae at visualized as concentric circles. 
2 : : This observation is explored 
(There is no need to include arbitrary constants when integrating 2y and —2x Pert semen at 


since they have already been taken care of by the constant C.) Since this last programme. 


22 MST204 2.3 


equation shows how the variables x and y are related for each value of C, it 
implicitly describes the family of solutions of our differential equation. However, 
solutions of differential equations are functions and so we must find a formula 
expressing y explicitly as a function of x. This can be done by treating the above 
implicit solution as an algebraic equation and solving it for y to obtain 


y= +./C — x’. 
(The + sign means that for each value of the constant C there are two distinct 
functions satisfying the differential equation: y = op C — x? and y = —./C — x’.) 
Notice that in this solution the constant of integration C is no longer added to the 
right-hand side of the equation but has been absorbed under the square root sign. 
Nevertheless, by giving C its various values we still generate the whole family of 
solutions. The above formula is therefore the general solution of the differential 
equation. 


Since the road by which we arrive at such solutions is a fairly long one, the results 
should always be checked by substituting in the original differential equation. In 
the present case, differentiation gives 


dy ae x Wherever the symbols + or 
dx ee + appear, the upper sign 
Jc . refers to the solution 
so that y =./C — x’ and the lower 


sign to the solution 


d 
ay = (+2,/C — x?) x [ss y= —./C-x’ 
—x 


which agrees with our original differential equation. 


The procedure described above can be summarized as follows. 


Procedure 3.2 
1. This procedure applies to differential equations which can be 
written in the form 


0) =f) 
X 


where f and g are given functions. 


2. Integrate both sides with respect to the same variable x 
obtaining 


d 
Jew) geae Ee [reo dx + C. 


Always be sure to include the arbitrary constant C. 


3. The rule for integration by substitution justifies writing this as 


[eonay = [reyax +c. 
4. Carry out the two integrations, to obtain a formula relating x 
and y. (This assumes, of course, that the integrations can be done.) 
5. Rearrange the formula in Step 4 to give y in terms of x. 


6. Check that the resulting formula for y satisfies the differential 
equation. If it does then you have found the general solution of 
your differential equation. 


Exercise 2 
In each of the following, ‘solve’ means ‘find the general solution of’. In part (iv) the 


1 
condition (y > 0) ensures that — makes sense; it means that you are only required to find 


solutions for which y is positive for all values of x. 
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(i) Solve (111) Solve 
d 
yo =e" Mees x 
dx 1 + y* dx 
and check your solution. and check your solution. 
(i1) Solve (iv) Solve 
dy 1 1 dy 
[eae —-— =] >0 
dx 1+%x? y dx VA 
and check your solution. and check your solution. 


[Solution on p. 45 | 


3.3 Separation of variables 


This subsection is about a type of differential equation which arises very frequently 
and which can be solved by bringing it to the form considered in the previous 
subsection. One equation of this type is 


a a (1) 
dx 


This cannot be integrated directly because the variables x and y appear together 
on the right-hand side, however, we can bring it to the form 


dy 
B(y) 7 = F(*), 
ik 
which can be integrated using Procedure 3.2, by dividing both sides of Equation (1) 
by 1 + y*. This gives 


Le 
1 + y* dx 


which is one of the equations you solved in Exercise 2, obtaining the general 
solution 


y = tan (4x? + C). 


This must also be the general solution of Equation (1) for we do not alter the 
information contained in Equation (1) by dividing it on both sides by 1 + y?. 


The procedure we have used here is known as separation of variables, because by 
dividing by 1 + y* we brought the equation to a form where the left-hand side 
depended only on the variable y and its derivative dy/dx, while the right-hand side 
depended only on the variable x: the two variables x and y were separated to the 
two sides of the equation. The procedure can be summarized as follows. 


Procedure 3.3: Separation of variables 


1. This procedure applies to differential equations which can be 
brought to the form 


d 
ae = SOM(Y) 


where f and h are known functions. 
2. Divide both sides by h(y), obtaining 


1 dy 
h(y) dx = f (x). 


3. Apply Procedure 3.2 to this last equation. 


A complication which sometimes arises is that h(y) may take the value zero for 
some ‘critical’ value or values of y. Since division by zero is meaningless, a solution 
obtained by the separation of variables method can only be relied on to the extent 
that it avojds such critical values of y. 
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If there are critical values of y (that is, values of y for which h(y) = 0) then there 
will be solutions of the differential equation which cannot be found using the 
separation of variables method. For example, any constant function, with constant 


os dy 

equal to a critical value of y, will satisfy the equation by making both = and h(y) A constant function is one of 
Z€TO, but to find such a solution using separation of variables would involve = a aa Paget - . 
dividing by h(y) = 0. However, there will also be solutions which avoid the critical constant value C for all values 
values of y and they can be found using the separation of variables method, as the of x. 
following example illustrates. 
Example 1 
Consider the equation 

dy 

a a 2 

Fa (2) 


This has the standard form of Procedure 3.3, with f(x) = 1 and h(y) = y. In this 
case the only ‘critical value’ of y is zero. One solution of the differential equation, 
not given by the separation of variables method, is therefore the constant function 
y = 0 (that is, the function which takes the value 0 for all values of x). All the 
other solutions can be found using the separation of variables method, but to 
avoid the critical value 0 we treat solutions for which y > 0 and y < 0 separately. 


Dividing both sides of the equation by y (valid for all non-zero values of y) gives 


so that 
1 
{ -dy=x+C (3) 
y 
where C is an arbitrary constant. Now we want a primitive of 1/y. If y > 0 we can 
use the obvious one, which is log, y, obtaining 
log.y=x+C (y > 0). 
Solving for y then gives 
y an etc. 


But if y < 0 the primitive log, y will not do since negative numbers do not have 
(real) logarithms. The appropriate primitive of 1/y for negative y is log, (—y). (If 
you have not come across this before you should check that it does have 1/y as its 
derivative, and also note that —y is positive so that it does have a logarithm.) 
Using this primitive in Equation (3) gives 


log. (—y)=x+C (vy < 0), 
and solving for y then gives 
y= —ertc 


So it seems that the family of solutions of the differential equation requires the 
three formulae 


yee, y= ey = (4) 


for its complete description. However, in Section 4 we shall solve the same 


differential equation by a different method and obtain the general solution as a 
single formula: 


y= Cle (5) 


where C is an arbitrary constant. On the face of it this formula looks quite 
different from Formulae (4) but, as we shall see in Section 4, Formulae (4) and 
Formula (5) describe the same family of solutions. We shall return to this in the 
tape commentary for Section 4 but in the meantime you may like to think about it 
yourself. 
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Exercise 3 


Obtain and check the general solutions of the following differential equations. (In some 
cases the allowed values of x and y have been restricted to ensure that you will not be in 
danger of dividing by zero.) 


dy 

(i) = 

= dy 

(i) x—=y (x > 0,y > 0) 

dx : 

d 

Gi) 2=- (y>0) 
i 


(Solution on p. 45] 


Exercise 4 


Obtain and check the general solution of the following differential equations, which are 
used in a later unit. In each case A and B are positive constants. 


dy : dy 
(i) ge (y > 0) (iv) A--=B-—y (B-—y<0O) 
x dx 
——— dy 
Sees ae PS age 2 
(11) re y (y > 0) (v) A = B’+y 
dy 
(ii) A—=B-y (B—y>0) 
dx 7 


[Solution on p. 46] 


3.4 Partial fractions 


In Unit 3 we shall need the solution of the following differential equation 


dy . 
Pi ea 


where a and b are constant. Provided that ay — by* # 0 (that is, y 4 0, y a/b) we 
can use the separation of variables procedure to bring the equation to the form 
1 dy 
ay = by dx 


and integration on both sides then gives 


dy 
Se ek ae C 
{- ei bp xX + 


where C is an arbitrary constant. But to complete the solution we need to evaluate 
the integral on the left-hand side, which is of a type you may not have seen before. 


As always with integrals, we try to transform the integral to one of the types we 
know how to deal with. In this case neither integration by parts nor substitution 
(unless you are extremely clever) is any help; but the following procedure known 
as a partial fraction expansion, solves the problem. 

We notice that the integrand has a denominator which can be factorized: 


1 = 1 
ay—by* (a—by)y 
The trick is to look for two numbers N, and N, with the property that 


1 N N 
——_—— =—+_+-# (y¥#0Oand y ¢a/b). (1) 
(a~Dy)y @—ey -¥ 
Of course it is not obvious a priori that such numbers exist; but if they do then we 
can multiply both sides of the last equation by (a — by)y to obtain the equivalent 
condition 


1=N,y+N,(a — by) 


(y #0 and y a/b) 
that is 


1=N,a+(N, —N,b)y (y #0 and y £ a/b). 


Hint: et” = &e’ 
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This condition is satisfied for all values of y (except 0 and a/b) if and only if 
1=N,a and O0=N, —N,b. 

This last condition can be solved for N, and N, to give 
N,=1/a and N, =N,b=D/a. 


So provided that a is different from zero our supposition that the numbers N, and 
N, exist is justified and we have, from Equation (1) 


1 So l/a 
(a~—yyy aby 
b 


1 
= + (y £0 and y £a/b). 
a(a—by) ay 


It is a good idea to check such equations by trying a particular value of y; for 
example, substituting y = 1 into the right-hand side gives 
b 1 b+(a—b) 
a(a — b) = a(a — b) 
1 
ae: 


which is the same as the expression obtained by substituting y = 1 into the left- 
hand side. 


Now we can do the integration. To avoid the critical values 0 and a/b it is 
necessary to consider various ranges of y separately. We start by considering the 
case y>0,a—by>0: 


J" 5 = [E+ 
ay—by? J \a(a—by) ay ‘ 


b d 1 (d 
=-| Yt. 
a—by aJy 


1 1 
= —-log, (a — by) + —log.y (since y > 0, a — by > 0) 


and so we can finish solving our differential equation. We had already got as far 
as 


dy 
for x + 


we can now write 


1 1 
——log, (a — by) + -log, y=x+C. 
a a 


To complete the solution we want to solve for y. Multiplying both sides by a and 
using the properties of logarithms gives 


y 
ie a 
08.7 by a(x + C) 


taking the exponential of both sides and multiplying by a — by gives 
y = (a— by)exp [a(x + C)] 

hence 
y(1 + bexp [a(x + C)]) = aexp [a(x + C)]. 

Finally we have 


_ aexp [a(x + C)] 
~ 14+ bexp [a(x + C)] 


a solution which is valid provided y > 0 and a — by > QO, that is, provided 
0 < y < a/b. Other ranges of y can be treated by the same method. 
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A similar method can be used for any integral where the integrand is the ratio of 
two polynomials. The procedure described below covers only the simplest cases to 
which the method applies. One example of a more complicated case (with a cubic 
in the denominator) is given in Exercise 5(iv). 


Procedure 3.4: Partial fractions 


To express in partial fraction form the function 


ax + B 
as ie +c 


where a, b, c, «, 6 are numbers (with a ¥ 0) such that the 
denominator can be factorized into distinct factors. This condition 
will be satisfied if b* — 4ac > 0; for then the quadratic equation 
ax? + bx + c = 0 will have two distinct real roots. 


1. Let x,, x, be the solutions of ax? + bx + c = 0. Then the 
function becomes 


ax + B 
a(x — x,)(« — x2) 


2. Suppose that numbers N,, N, exist such that 


ax + B i - N, 
a(x —x,)(x-—x,) x-—x, 5 eee 


(x #x, and x #x,). 


3. Multiply both sides by a(x — x,)(x — x,) and collect terms: 


ax + B = aN, (x — x.) + aN,(x — x,) 
=a(N, + N,)x — a(N,x, + N,x,) 


4. Since this last equation holds for all x (except x, and x,) we 
must have 


a=aN,+aN, 


B = —ax,N, — ax,N . 


5. Solve the above pair of simultaneous equations for N , and N,. 


6. Check that your result 


dee Sees = N, 
ax + bx +e X=X, X= Xe 


is correct by substituting one or two particular values of x (not x : 
OF X51. 


This procedure does not cover the cases where the quadratic has one or no real 
roots. If it has one real root, say x,, then the substitution x — x, = w will simplify 
the integral. If it has no real roots then a substitution involving arctan will 
simplify the integral (cf. solution to Exercise 4(v)). But such integrals are beyond 
the scope of this section. 


Exercise 5 


Find partial fraction expansions for 

1 

Sas 
x+1 


(u) s +57 


1 ea 
ill) (where B is a positive constant) 
B2 es y? 


1 
aie mae 
[Solution on p. 47 | 


Exercise 6 

Using the partial fraction expansions you 
equations 

eee 

ie Pare Bae (y > 0) 

ee x+1 

(11) i See (-—2<x <1) 
Ss dy ae 

(i11) oF —y (-—B<y<B) 


[Solution on p. 47] 
Summary of Section 3 


Differential equations _ 
The solution of the equation 


d 


oe 
Pes 


with f a given function, is 
y= { S(x)dx +C 


where C is an arbitrary constant. 


To solve the equation 
d 

(v= =f) 
x 


make y the subject of the formula 


few dy = {r (x) dx + C. 


To solve the equation 
d 
= = SHY) —(h() #0) 
x 


make y the subject of the formula 


1 


The above are general solutions involving a constant C; to obtain the particular 
solution satisfying a given further condition, use the further condition to evaluate C. 


Integrals 
| dx = 
aed 


To integrate the function 


1 
—log, (ax + b) 
a 


ax + B 
ax? + bx +c 


‘a as — b) 
a 
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found in Exercise 5, solve the differential 


ifax +b>0 


ifax+b<0O 


where b* > 4ac, first express it in partial fractions: 


i Ne 


where N,, N,, are constants and x,, x, are the roots of the quadratic ax* + bx +c. 


You will meet this equation 
again in Unit 4 


Procedure 3.2. 


Procedure 3.3. 


Procedure 3.4. | 
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4 Linear equations 


4.0 Introduction 


There is no universal method for obtaining the general solution of a differential 
equation, but the method to be described in this section, the integrating factor 
method, does solve an important class of equations, called linear equations, for 
which the separation of variables method may not work. Like the separation of 
variables method, the integrating factor method depends on rearranging the 
differential equation to a form that can be integrated directly. We begin by 
discussing the type of direct integration appropriate for this. 


4.1 Direct integration using the product rule 


The form of direct integration considered here is based on the formula for 
differentiating a product. The product we are interested in is p(x)y, where p is 
some function to be chosen later. The formula for differentiating p(x)y is 


ata 
dx 


d = dy 
qn) = 00 


Consequently, if the differential equation we are trying to solve can be written in 
the form 
dy _ dp(x) 


ies dx 


p(x) y = q(x), (1) 


for some functions p and q, then the differential equation is equivalent to 
d 
fe ae) 
x 


and can therefore be integrated directly, to give 


p(x) = | olds £6 


Dividing through by p(x) then gives y as a function of x. 


So any differential equation that can be brought to the form of Equation (1) can 
be integrated directly. The next subsection describes how this can be done for 
certain equations. 


Exercise 1 


Check that the following differential equations have the form of Equation (1). Determine p 
and q in each case, and hence integrate the equations. 


: dy 
ae a4 
(1) par eae 


Sa 
(ii) sinx— + (cosx)y = sinx 
dx 


d 
(iii) en + 3e*y = e* 
x 


[Solution on p. 48 | 


4.2 The integrating factor method 


The method we are going to study here applies to any differential equation of the 
form 
d 
a = A(x) + Iex)y (2) 
x 
where k and / are given functions. An equation having this form is said to be 
linear. 


There is an analogy here with 
linear first order recurrence 
relations, that is, recurrence 
relations of the form 


Ung, =k, + Lu,. 
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To solve Equation (2) we want to find a way of rearranging it to the form 


ay og 
poy + Py = ax (1 


which we saw how to integrate in the previous subsection. (We do not yet know 
what the functions p and q are, but we shall find out during the course of the 
calculation.) As a first step towards bringing Equation (2) to the form of Equation 
(1), we subtract [(x)y from both sides of Equation (2), obtaining 


d 
se — Holy = k(x) (3) 


To make the first term of this equation match the first term of Equation (1) we 
multiply both sides by p(x), to get 


d 
p(x) — pOxllely = plx)k(x) (4) 


The remaining terms of this equation can be made to match Equation (1) by 
choosing the functions p and gq so that 


d 
Pe) = —¢xplee ) 


and 


q(x) = p(x)k(x). 


In fact these two conditions are sufficient to determine the two functions p and g. 
We are particularly interested in the function p, because multiplication by p(x) 
takes us from Equation (3), which is not of the directly integrable form, to 
Equation (1), which is. For this reason p(x) is called an integrating factor. 


The condition determining p(x) is Equation (5), a differential equation which can 
be solved by the method of separation of variables. Dividing both sides by p(x) 
and then integrating both sides gives 


1 dp(x) 
p(x) as dx = | (-t09) dx 
that is 
dp(x) 
[ne = | (-teoyae 


and so one solution is given by 


log, p(x) = | (=IGx)) dx 


hence 


p(x) = exp) | (Ibo) de | (6) 


No arbitrary constant is necessary, nor need we look for solutions with p(x) < 0, 
because we only need one way of bringing our differential equation to the form of 
Equation (1), not all the possible ways. Having determined the integrating factor 
p(x) from Equation (6), we can solve Equation (4) (and hence our original 
differential equation) by direct integration, after which it is easy to find y. The 
following example should clarify the procedure. 


Example 1 

Let us apply the integrating factor method to the equation 
d 
=x 4iy. (7) 
dx 


This has the linear form of Equation (2), with k(x) = x and I(x) = 1. The first step 
is to subtract the term containing y, which in this case is y itself, from both sides, 
obtaining 


Equation (7) is discussed in 
the television programme for 
this unit 
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Soy wall 


Now we want the integrating factor p(x). This is best found from Equation (6): 


p(x) = exp [(-te ix| 
= exp} f(—tyae | 


x 


=e”, 
Multiplying both sides of Equation (8) by this integrating factor we obtain 


(== =f ye x. (9) 


If the integrating factor has been calculated correctly, the left-hand side of 
Equation (9) should be the derivative of the product p(x)y, that is of e *y. Using 
the product rule to calculate the derivative of e~*y we find that Equation (9) can 
indeed be written 


ges 
= 1e “Pi = ¢- * x. 
cree SE 
Direct integration now gives 
¢ y= fers dx 
= —xe *+ | e dx (integrating by parts) 


= —xe *~*—e *+C. 


Finally we solve for y, to obtain the general solution of our Differential Equation The solution given in the 
(7): introduction to this unit is a 
particular case of this, with 
y= Ce —-x-1. E> 1. 


The general procedure for solving a first order differential equation by the 
integrating factor method is given as Procedure 4.2 on the next page. 
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Procedure 4.2: The integrating factor method 


1. This procedure applies to differential equations which can be 
written in the form 


dy | 

oe k(x) + I(x)y 
where k and / are given functions; such equations are said to be 
linear. 


2. Subtract /(x)y from both sides, obtaining 


dy = 
Bes I(x)y = k(x). 


3. Define a function p (the integrating factor) by 


p(x) = exp] (—I(x)) éx| 


(no arbitrary constant is necessary). 


4. Multiply both sides of the equation in Step 2 by p(x), 
obtaining 


plo) — plxplbe)y = pla k(x) 
= 4 


5. Verify that the left-hand side is the derivative of p(x)y (if not, 
you have made a mistake in calculating the integrating factor) 
and integrate both sides, obtaining 


p(x)y = [r(K(x) dx + Constant. 


Do not forget the arbitrary constant. 
6. Carry out the integration of p(x)k(x). 
7. Solve algebraically for y (that is, divide both sides by p(x)). 


8. Check that the expression you have obtained for y satisfies 
your original differential equation. If it does, you have found the 
general solution of your differential equation. 


4.3 Using the integrating factor method (Tape Subsection) 


The purpose of this subsection is to give you practice in using the integrating 
factor method. It consists of three examples accompanied by a tape commentary. 


d 
The first example, the equation = = y, has already been solved in Subsection 3.3 
2 


using the separation of variables methed; but the method used here is better 
because it covers the cases y > 0, y = 0 and y < 0 all together instead of 
separately. 


While working through the examples you will need to refer to the steps in 
Procedure 4.2. 


Start the tape when you are ready. 
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¢ A worked excamplie 


a The integrating factor 


exp [ ((-Ucxy do | 


exp [ | (- eae 


& Rearrang! ng 


& Solve and check 


Step 7 4 


d 
Gri & 
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© Less help 


P(xc) = excp | S- ( ) dx | 


C ) 


= 
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ee ee en, aa - Y= = 
Ax 


Manipulate 
+ identify 
Steps 1,2 


Manipulate 
steps 4,5 


In cegrate 
Step 6 


Solution y = x log =< +Ax 
Step 7 


© Solution? to dy =y 


36 MST204 2.5 


d 
Solution to oa y 
dx 


At the end of the tape commentary we asked you to reconcile the general solution 
y = A& (10) 
obtained by the integrating factor method, with the general solution 
y=ete, pees y=0 (11) 


obtained using separation of variables. To do this observe that positive values of 
the arbitrary constant A in Solution (10) can equally well be written e“ to give the 
first of Formulae (11), whereas negative values of A can be written as —e“ to give 
the second of Formulae (11). Setting A = 0 gives the remaining solution y = 0. 


The above example illustrates the fact that it is often possible to find a more 
economical way of writing the general solution of a differential equation simply by 
redefining the arbitrary constant. 


Summary of Secticn 4 
To solve the linear equation 


OY = k(x) + Mx) 


XxX 


subtract I(x)y from both sides and then multiply both sides by the integrating 


factor 
p(x) = exp| [- lode 


The resulting equation has the form 


dy ad 
ply + Py = a(x) 


and integrates directly to give the general solution 


p(x)y = { g(x)dx + C. 


5 Revision 
5.1 Choosing the best method 


This unit contains a number of different methods for solving differential equations. 
You will want to know how one decides which to use. 


The first step is to bring the equation to the standard form 


d 
= — m(x, y). 


Not every differential equation can be brought to this form, for example the 
equation may include second or higher derivatives, but such equations are beyond 
the scope of this unit. 


Next look to see if an exact solution is possible by one of the methods described 
in Sections 3 and 4. In effect, these methods cover two cases only: 


(i) equations soluble by the separation of variables method (Procedure 3.3), 
d 


yt 


(ii) linear equations using the integrating factor method (Procedure 4.2), 
d 
= k(x) + U(x)y, 
dx 


where f and h, or k and |, are known functions. The various cases of direct 
integration can be regarded as special cases of (i) and (ii) (for example, setting 
h(y)'= 1 or I(x) = 0 gives the direct integration considered at the beginning of 
Section 3). 


Equations containing second 
order derivatives are discussed 
in the next unit on differential 
equations. 
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If the exact methods fail, either because m(x, y) has the wrong form or because you 
cannot do the integrations, or even because (as can easily happen) the formula 
given by the exact solution method is too complicated to be useful, you will have 
to fall back on an approximate method. This usually means a numerical method 
and for the time being, therefore, the Euler method. By taking a fairly large step 
size you can carry out the Euler method using a calculator, but if you want 
anything better than a qualitative picture of the solution, you will have to use the 
computer. 


A disadvantage of the Euler method (or any other numerical method) is that it 
only gives particular solutions. If you require information about the general 
solution for a differential equation you cannot solve exactly, you can obtain it by 
plotting the direction field and then sketching enough solution curves to give a 
good picture of their typical behaviour. 


Exercise 1 


For each of the following equations, select the option(s) which give(s) the most appropriate 
method(s) of solution. 


A direct integration 
B separation of variables 
C the integrating factor method 


D none of these 


dy | dy : 
(1) ees (vi) oe 
s2 dy x dy . 
(11) = y (vi1) a x*+y 
(iii) — (viii) oe 
dx ! dx . 
(iv) ees (ix) o 
dx y dx 
dy r 
(v) eis y 


[Solution on p. 48] 
5.2 Further exercises 


The exercises in this subsection will enable you to revise and consolidate what you 
have learned from the rest of the unit. 


Exercise 2 
For the direction field shown, sketch the trajectory through the point P. 


[Solution on.p. 48] 


Exercise 3 

For the differential equation 
= =1-—-xy 
dx 

with y(0) = 1, 


(i) write down the recurrence relation for Euler’s method with step length 0.1; 


(ii) using step length 0.1 calculate the Euler’s method approximation to y(0.3). 
[Solution on p. 48 | 
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Exercise 4 
If the step size in the above calculation were reduced to 0.05, 


(i) would you expect the error in y(0.3) to increase or decrease, and by what factor 
(roughly)? 


(11) would the labour of calculation increase or decrease and by about what factor? 
[Solution on p. 48 | 


Exercise 5 (Do not attempt this exercise before studying Section 6.) 


If the step size were further reduced to 0.001, write down the options and responses you 
would use to carry out the calculation using the computer program RECREL with print- 
out of the approximations to y(0.1), y(0.2) and y(0.3) only. 


[Solution on p. 48] 

Exercise 6 

The general solution of a certain differential equation is 
y= Ox 4. 

Find the particular solution for which y(4) = 8. 

[Solution on p. 48 | 


Exercise 7 
Find general solutions of the differential equations 


2 
oe oe nee 
PS ae 
dy 2 
Gy =o 225 “228 
qe = 


(i) =y+y? >) 


(iv) = 2y—x 
[Solution on p. 48 | 


Exercise 8 
Carry out the integration 


x+4 
—_____—_—d —2). 
lee oe 


[Solution on p. 49 | 


6 Computer terminal visit for Units 1 and 2 


6.0 Introduction 


See the course guide and possibly also the Stop Press for up-to-date information 
about this visit. Be sure to check whether there is a TMA question requiring use 
of the terminal. 


The purpose of this visit to the computer terminal is twofold; (i) to enable you to 
learn how to use the computer package RECREL to make rapid calculations 
involving recurrence relations and Euler’s method, and (ii) to see some of the 
numerical phenomena described in Units 1 and 2, such as ill-conditioning and the 
balance between accuracy and cost, in an actual computation. 


Before going to the terminal you should have prepared yourself by studying 
Section 6 of Unit 1 and this section. Section 6 of Unit 1 described in general terms 
how to work the computer packages, and showed you how to turn algebraic 
expressions into valid input expressions for the computer. This section contains a 
description of the package RECREL and some home exercises to help you 
familiarize yourself with it. After completing the home exercises you will be ready 
to start planning your work at the terminal based on the computer exercises in 
‘Subsection 6.3. It is essential that you prepare these computer exercises carefully at 
home so that you will not waste time at the terminal. 
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6.1 The computer package RECREL 


The computer program RECREL is a package to use when solving first and second 
order recurrence relations. The package can also be used to implement Euler’s 
method to solve first order differential equations. An example of the use of this 
package was given on p. 44 of Unit 1. You may need to refer to this to help you 
assimilate the package. : 


Command options 


OPTIONS this tells the program to print a list of available options. 

SOLVE this tells the program to run the problem after all the data has 
been input. If some data is missing an error message will be 
printed. 

HELP this tells the program that you are stuck and need advice. 

LIST this tells the program to print the problem and data previously 
input so that you can check which problem you are solving. 

STOP this is the only way to stop the program. 

ROOTS for constant coefficient second order recurrence relations this 


option tells the program to compute and print the two roots of 
the auxiliary equation. 


ANSWER this gives the solution to the computing exercises so that you can 
check that your answer is correct. To the question 


EXERCISE? 


you respond with the number of the exercise you are interested in. 


Problem options 
10 Order of the recurrence relation 
To the question 
ORDER =? 
you respond with either 1 or 2 depending on the order of the recurrence relation 
that you wish to examine. This option is not essential for Euler’s method. 


11. Input recurrence relation 
To the question 


U(R +1) =? 


you respond by typing in the right hand side as a valid input expression (see 
Subsection 6.3 of Unit 1). For example, if you wish to examine the recurrence 
relation 


u,., = 3u, — 2u,_, + 3r 
you would respond by typing 

3*U(R) — 2*U(R — 1) + 3«R. 
12 Input m(x, y) for Euler’s method 
To the question 

Y=? 


you respond with a valid input expression for m(x, y). For example if the 
differential equation were 


dy 
cs ee 3 
= xy + 


your reply to Y’ = ? would be 
XY + 3. 


y’ 1s another way of 


writing = 
a 
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Method options 

20 Forward recurrence 

This method uses forward recurrence on either a first or second order recurrence 
relation. The number of terms in the sequence is specified using Option 35. 

21 Backward recurrence 


This method uses backward recurrence on a linear first order recurrence relation 
using Option 32 to specify the starting point. The program will automatically 
adapt the recurrence relation specified by Option 11. For example if you want to 
examine the recurrence relation 


Ce 1 ae es 1)u, 


then you input this recurrence relation using Option 11. The package will 
automatically convert this into the form 


l—u 
u, = ——_+ 
r+1 


for use with backward recurrence. 
22 KEuler’s method 
This uses Euler’s method to solve the differential equation 


= = m(x,y) 


where m(x, y) is specified using Option 12. 


Parameter options 

30 Initial condition u, 

This parameter specifies the initial condition u, for first and second order 
recurrence relations. 

31 Second initial condition u, 

This parameter specifies the second initial condition for second order recurrence 
relations. 

32 uy for backward recurrence 

This is used to input the initial condition for backward recurrence with a linear 
first order recurrence relation. N is specified by Option 35. 

33 x, for Euler’s method 

This is used to input the initial value of x for Euler’s method. 


34 yo for Euler’s method 
This is used to input the initial value of y for Euler’s method. 


35 The number of terms /N in the sequence 

This is used to input the number of terms in the sequence. This must be specified 
for both forward and backward recurrence and for Euler’s method. 

36 Step size h for Euler’s method 

For Euler’s method this option specifies the step size. 


Print options 

40 Outline printout 

If you only want selected output you should use this option to specify k such that 
only every kth term in the sequence is printed. 

41 Full printout (default) 

All terms in the sequence will be printed. 


A default option is the one the 
program will use if no other is 
specified 
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42 Print solution only 


Only the last term in the sequence is printed (intermediate terms will not be 
printed). For backward recurrence the value of u, will be printed. 


6.2 Home exercises 


The following exercises are to enable you to familiarize yourself with RECREL so 
that you will not waste time at the terminal. 


Exercise 1 


(i) It is proposed to put the calculation described in Subsection 2.2 on to the computer; 
that is, to calculate y(1) given that 


dy 
Pie oS 


and y(0) = 0, with h = 0.2. Write in the table the options you would use and the responses 
to each (where called for) in order to input this problem and print the calculated values of 


y(0.2), y(0.4),..., y(1). 


Information Response 


Formula for m(x, y) 
Method 

Intitial value of y 
Initial value of x 
Number of steps 
Step size 

Pes of printout 


(11) Suppose you wanted to increase the accuracy by a factor of 100, so that h is reduced to 
0.002; but that to save printing time you only wanted to print approximations to y(0), 
y(0.1), y(0.2),..., y(1). What changes to the answer to (i) would this entail? 


[Solution on p. 50] 


Exercise 2 
In Sections 3 and 4-of Unit 1 we considered the recurrence relation 
uns — Nt ares (r 5 1)u, 


with uy, = 1 — e | = 0.63212056..., which can be used to calculate a sequence of values for 
the integral 


1 
‘= } ot ae: 
1) 


Make a table, like the one in Exercise 1, showing the information you must give the 
computer, the options you would use, and the responses to those options, if you wanted to 
calculate the first 25 terms of the sequence, using 


(1) forward recurrence 
(ii) backward recurrence, starting at u,. 


[Solution on p. 50] 


6.3 Computer exercises 


Now here are the exercises for you to try on the computer. Before setting out, 
decide on the order in which you will try them; then decide, for each problem, 
which options you will use and which responses. In this way you will make the 
best use of your very limited time at the terminal. Solutions to these exercises can 
be found from the package using the command option ANSWER 


Note: There may also be a TMA question to be done at the terminal. 


Exercises 3—6 repeat the same numbered exercises from Subsection 6.4 in Unit 1, and 
are included here to refresh your memory. Exercises 7 and 8 you have not met before. 


42 MST204 2.6 


Exercise 3 
Use the package to carry out the problem described in Exercise 2 (use both forward and 
backward recurrence). 
Exercise 4 
The recurrence relation 
ta a 
with uy = e — 1 = 1.7182818... can be used to evaluate a sequence of values for 
1 
i= | e"* be 
0 
Use the package to compute the first 25 terms of this sequence. 


Exercise 5 


Check that your mortgage repayments are correct so that you repay the loan in the 
specified time period. (Note that your mortgage may include property and/or life insurance 
which are extras and should be omitted from the calculations.) 


If you do not have a mortgage, work out the monthly payment for a £30,000 mortgage 
borrowed over 25 years if the interest rate is 17°. Use the package to determine the 
amount outstanding at the end of each year. 


Exercise 6 


Use the package to determine the behaviour of the sequences generated by the following 
second order recurrence relations. 


(1) Up+1 = Uy — Uy-1, Up = 0, uy = 1 

(it) 444 =U, — 2u,-1, i = 0,0; = 1 

(ii) up4, = 0.9u, — 0.2u,_, + 1500, Up = 3600, u, = 4000 
(iv) u,+1 = 8u, — 16u,_, + 4, Up = 3, U, = $. 


Exercise 7 
Use the package to carry out the problems described in Exercise 1 (i) and (ii). 


Exercise 8 


This exercise suggests a possible way of treating problems involving differential equations 
which you can not solve exactly. 


(i) Obtain a numerical solution for the reflector equation The reflector problem was 
treated in the television 
ete te programme. 
x x 


which satisfies y(1) = 0, using a convenient step size, say h = 0.2. Take your solution as far 
as x = 3. 


(ii) Plot y against x forl <x <5. 


(iii) Try some other step sizes, including h = 0.1, h = 0.05, h = 0.02, and plot a graph to 
show how the computed value of y(5) varies with h. Estimate the exact value of y(5) from 
this graph. 


(iv) (Optional) What would happen if you used a negative h, with the same initial 
condition y(1) = 0? 
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Appendix: Solutions to the exercises 


Solutions to the exercises in Section 1 


a CG ou ee G8 
i, eer 3 
x x 
x? = x* Ear 3 
i ee 
x 


(c) 


4. 
the line segments. 
2. In each case we start by differentiating both sides of the 
given formula 
a -@ 
i |: ae: = 
(1) Ax = a & + 5) =3(2x) + 0=-x. 
a ae 
(ii) a qx 2 + C)=3(2x)+0=-x. 
(111) We have SGkUMEES BUSH ECEee 
ed = Ses —x—1)=Ce-1 From the diagram, the required family of curves appear to be 
dx dx straight lines with slope 2, so the family can be described by 
alten ety ox 4 (Ce — x — 1D) = Ce 1 the equation 
SO. —— st + y. 
dx 


d 
Differentiation gives = = 2+ 0 =2, which checks. 
x 


5. The function giving y on the lower half of the circle is 


y = —. /R? — x’. 


3. (i) Using the suggested systematic procedure we have 
Hence the slope of the circle at (x, y) is 


ay a = 
ie 2)1/2 _ 1(R2 _ 2 1/2 — 2x 
(a) = = —_(R? — x)? = A(R? — x?) 1? x (—2x) ee ae 
aoe dx 
er ee eS 1 
Ee ee ees 
x —Xx 2y 
(b) y = (R?2 — x?)1/?2 poe 
6. In each case the differential equation is found by 
d 
(c) so a = = equating ce with the slope of the direction field at (x, y) 
dx y dx 
ees 
(ii) The procedure used in part (i) gives ae 
iy & x wt OY 
(a) ee (ii) ss 


the line OP has slope and so the direction field at P also 


has slope = Hence the differential equation is 
x 


dy _y 


dx x 


MST204 2 Solutions 


Solutions to the exercises in Section 2 
1. 


+ 


a 


~~ i i ae 


~® CR SA TE Se 
Po ee a 


0.2 04 0.6 0.8 
— OO 


Ps 
5 ee 
~ CK OSA 


~~ ~~ 


The first line segment has slope 0 and ends at the point P, 

= (0.2,0). The next is drawn with a slope equal to that of the 
direction field at P, (the slope can be estimated from the 
diagram of the direction field, there is no need to calculate it) 
and continues until it reaches the point P, at which x = 0.4 
(i.e. until it meets the vertical line x = 0.4). The next has 
slope equal to that of the direction field at P, and continues 
until it meets the vertical line x = 0.6, and so on. 


2. The complete table is: 


Vr+i = yr x hm, 


0.04 
0.128 
0.2736 
0.48832 


The co-ordinates of P; are shown in bold type. 


3. The only change to the column headings is that the 
formula for m, is now y, instead of x, + y,. Also, since the 
initial value y(0) is now 1 instead of 0, the first entry in the y, 
column is now 1. The table is therefore: 


Vr+1 ee yr ¥ hm, 


0.41472 | 2.48832 


The approximation to y(1) is shown in bold type. 


4. (i) The total change in x is 2 — 1 = 1, so with n steps the 
step size is 1/n. 


(ii) Euler’s formula for the differential equation can be 
written 


W(X,41) = Y(x,) + hx, 

so we have 

(a) ifn=1 then h = 1 giving 
y2)~ ypi)t+tixl=24+1x1l=3. 
so for 1 step y(2) ~ 3, 
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(b) ifn =2 then h = giving 

y4) ~ y(1) +4x1=24+4x1=2 
y(2) ~ y(1z) +2 x 197423 +9 x 17 = 34 
so for 2 steps y(2) ~ 3.25, 

(c) ifn=4 then h = giving 
yiz)~24+4x1=24 

y(1z) = 24 +4 x 14 = 276 
y(1z) ~ 276 +4 x 13 = 258 

y(2) ~ 2g +4 x 1g = 35 

so for 4 steps y(2) ~ 3.375. 


(iii) The exact value of y(2) is 4x? + 3 with x = 2, that is 
Wes $s. 


For n = 1 the error is 3.5 — 3 = 0:5. 
For n = 2 the error is 3.5 — 3.25 = 0.25. 
For n = 4 the error is 3.5 — 3.375 = 0.125. 


(iv) It appears that the error is proportional to h, in fact 
the error seems to be equal to Sh. 


(v) For an accuracy of 10~* we would need 4h = 10~¢, that 
is h = 2 x 10~*. By (i) the number of steps would therefore 
be 1/h = 5000. 


Solutions to the exercises in Section 3 


1. (i) Integrating the function 6x’ gives the general 
solution 


y=2x +C 


The required particular solution satisfies y(1) = 5, that is 
2 x 1° + C =5, which implies C = 5 — 2 = 3. The particular 
solution is therefore 


y = 2x? + 3. 


~2* gives the general solution 


(ii) Integrating the function e 
y= —se *+C. 


The required particular solution satisfies y(0) = 2, that is 
—te° + C = 2, which implies C = 2 + 4 = 24. The particular 
solution is therefore 


y = —te7-* + 24. 
(111) Integrating the function asin bx gives the general 
solution 


y= —=cos bx + C. 


The required particular solution satisfies y(0) = 0, that is 
~ 008 (0) + C = 0 which implies C = - The particular 


solution is therefore 


= —“cosbx +2 
ea ie 


2. We solve each equation by following the steps listed in 
Procedure 3.2. 


(1) Steps 1 to 3 tell us that we can integrate both sides of 
the equation to obtain 


[ra = [xta + C. 

Step 4: doing the integration gives 

yo = 1x5 4. C. 

Step 5: solving for y we obtain the general solution 


y= /3x°? + 3C. 
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Step 6: check: 
1 


vo 


4 


dy 
= = 3(2x° + 3C)~25 x 3x4 = 4— x 3x 


d 
hence ys = x*, 
dx 


(ii) Steps 1 to 3 tell us that we can integrate both sides of the 
equation to obtain 


1 
2 ae eee 
fe y \a pete 


Step 4: doing the integrations we obtain 

e’ = arctanx+C 

(arctan x means the angle between = and i whose tangent 
is x). 

Step 5: the general solution is therefore 


y = log, (arctan x + C). 
Step 6: check: 


d { 
= = 4(3x3 + 3C) "2/3 x 3x4 = 33 ox 
hence oe = eae 
ax 1+ x" 


(iii) Steps 1 to 3 tell us that we can integrate both sides of 
the equation to obtain 


1 


Step 4: doing the integrations we obtain 
arctan y = 4x? + C. 

Step 5: the general solution is therefore 
y = tan (4x? + C). 

Step 6: check: 

= = xsec*(4x? + C) 


1 ed sec? (4x? + C 
hence oo Sola 


1+y?dx 1+ tan? x? + C) 
=x (since sec? 9 = 1 + tan’ 6). 


(iv) Steps 1 to 3 tell us that we can integrate both sides of 
the equation to obtain 


1 
[-ay = ax +. 
Bo 


Step 4: doing the integrations we obtain 
log. y=x+C (since y > 0). 
Step 5: the general solution is therefore 


y =exp(x + C). 


Step 6: check: 
Lays 1 


= +C)=1. 
ydx exp(x+C) ee ) 


3. We solve each equation by first applying Procedure 3.3 
and then following the steps listed in Procedure 3.2. 
(i) After dividing by e” the equation becomes 
dy 
= 2 
2 dx 


which has the form described in Step 1 of Procedure 3.2: 
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Steps 2 and 3: integrating and using the substitution rule 
gives 


feva o fea + C. 


Step 4: doing the integrations we obtain 
—e%=e+C. 

Step 5: solving for y gives the general solution 

y= Ele OC) 

Step 6: check: 

(11) The equation can be brought to the right form for 
Procedure 3.3 by dividing by x giving 

es | 
de sik - x); 

and we can bring this to the form described in Step 1 of 
Procedure 3.2 by dividing by y to obtain 

idy 

a 


Steps 2 and 3: integrating and using the substitution rule 
gives 


1 1 
[-a = [-ax 
y x 


Step 4: doing the integrations we obtain 


log. y 
Step 5: solving for y gives the general solution 


= log. x +0 (since x > 0 and y > O). 
y = exp (log,x + C) 

= xe“. 
Step 6: check: 

8 noes 

‘ax 


(11) After multiplying by y (or dividing by 1/y) the equation 
becomes 


dy 
a 
which has the form described in Step 1 of Procedure 3.2. 


Steps 2 and 3: integrating and using the substitution rule 
gives 


[vay = [xax ee 


Step 4: doing the integrations we obtain 

dy? = 4x7 + C. 

Step 5: solving for y gives the general solution 
y=./x* +2C 


Step 6: check: ae 2Cy 1/2) 
p 6: chec Len = (x* 426) "ite + 26) 7) = 


4. (i) Dhtviding by y and integrating both sides gives 


[oa = [(-ayax e 


log. y= —Ax+C 
Solving for y gives the general solution 


y = exp(—Ax + C). 


hence (since y > 0). 


d 
Check: = —Aexp(—Ax+ C)= —- 
x 
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(ii) Dividing by y* and integrating both sides gives 


[ow = |(-ayax + €, 


hence — —Ax+C. 


Solving for y gives the general solution 


1 
AG 


j= 


Gua Se Soe 
ie ax Cy: ee 


(iii) Dividing by B — y and integrating both sides gives 


(e5°- Jase 


Doing the integrations (using the substitution u = B — y) we 
have 


A 
_ [Fau =x+C 

u 
sO —Alog,u=x+C (since u = B— y > 0) 
1.€. —Alog, (B-—y)=x+C. 


The general solution is obtained by solving for y: 
B —-x—C 
y p a 


hence y=B-—e *4e (4, 


d 
Check: Am = e~(/4e-*/4 = B ). 
dx 


(iv) The same steps as in (iii) brings us to 
_ [- du=x+C 
u 
where u = B — y. But since we now have u < 0 the integral 
gives 
—Alog,(—u)=x+C 
i.e. —Alog,(-—B+y)=x+C. 
The general solution is obtained by solving for y: 
—B+y=exp 4] 
A 


i.e. y=B+e *4e “4 


—C/Ag—x/A 


d 
Check: A = —e = B—y. 
dx 


(v) Dividing by B? + y’ and integrating both sides gives 


A 
i 5 dy = ax + C. 


Doing the integrations (using the substitution y = B tan u if 
necessary) we have 

A y 

—arctan— =x+C. 

B B 


Solving for y gives the general solution 


a eee 
:= rae 


d 
Check: yaa = B? sec? 
dx 


B(x + “ 
A 


= B? f + tan? PEO) 


— BR? a at 


MST204 2 Solutions 


5. We calculate each partial fraction following the steps 
listed in Procedure 3.4. 


(i) Step 1: factorize the denominator to give 
2 1 

eae ae 1 

Step 2: assume N, and N, exist such that 


x(x+1) x. x+1 


Step 3: multiply both sides by x(x + 1) and collect terms to 
obtain 


1=N,(x+1)+N,x 
=(N,+N,)x+N, 

Step 4: it follows that 

O0O=N,+N, 

ba, 


Step 5: this is a simultaneous equation with solution N, = 1, 
N, = —1, the partial fraction expansion is therefore 


1 1 1 


Step 6: as a check substitute x = 1 into this expansion; both 
sides then equal 4. 


(ii) Step 1: the solutions of x* +x —2=0Oare 1, —2. 
Hence 
x +1 x+1 
x 4x—2 (x—I(x+2) 
Step 2: assume N, and N, exist such that 
x1 N, N, 
(x aL ee 


Step 3: multiply both sides by (x — 1)(x + 2) and collect 
terms to obtain 


x+1=N,(%+2)+N,(x — 1) 
=(N, + N,)x + QN, — N,). 
Step 4: it follows that 
1lL=N,+N, 
1 =2N, —N,. 
Step 5: adding these equations together gives 3N, = 2, that is 


, = 4, substituting this into the first equation gives N, = 3, 
the partial fractions expansion is therefore 


x+1 aE 2 z 1 
XO Ae te 3 at 2) 


Step 6: as check substitute x = 0 into this expansion; both 
sides then equal — >. 


(iii) Step 1: factorize the denominator to give 
1 
a ee 
Step 2: assume N, and N, exist such that 
1 ay N, 

(B-y)(B+y) B-y Bry 
Step 3: multiply both sides by (B — y)(B + y) and collect 
terms to obtain 
1=N,(B+y)+N,(B-y) 

=(N, — N,)y + (N,B + NB). 
Step 4: it follows that 
0O=N,-N, 
1=N,B+N.B. 
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Step 5: adding B times the first equation to the second gives 
1 a = 
1 = 2N,B, that is N, = DB? substituting this into the first 


equation gives N, = aR? the partial fractions expansion is 
therefore 
1 1 1 


By BB) We +». 
Step 6: as a check substitute y = 0 into this expansion; both 
sides then equal 1/B?. 


(iv) The denominator of this expression is a cubic, but with 
minor modifications Procedure 3.4 can still be used. 
Factorize the denominator to give 


1 1 
se = xix — 1x 41) 
and assume N,, N, and N, exist such that 


1 eee N; 
= + + + : 
wie ee eb] 


x-1 


Multiply both sides by x(x — 1)(x + 1) and collect terms to 
obtain 


1 = N,(x* — 1) + N,(x? + x) + .N3(x? — x) 
=(N,+N,+N,)x?+(N,—N;3)x—N, 

Since this must hold for all x (except 0, 1 and —1) we have 

0O=N,+N,+N, 


1 = —N,. 
From the third equation, N, = —1, and from the second 


equation N, = N,. These two results, together with the first 
equation, imply that 0 = —1 + 2N,, that is N, = N; =}. 
The partial fraction expansion is therefore 


a 1 = 1 
ae 9) x eT ee eda 


As a check substitute x = 2 into this expansion; both sides 
then equal 2. 


6. (i) After dividing both sides of 
dy . 
oes Sas 


by y”? + y, we can integrate and use the substitution rule to 
obtain 


1 
ays fase 
or 


Using the partial fraction expansio1. »btained in 5(i), this can 
be written 


1 1 | 
can ree = {d oS 8 
IG il” : 


Doing the integrations, we obtain 
log. y —log.vVt1)=x+C 
By the properties of logarithms, this can be written 


(since y > 0). 


y 
log. ———— = x + C. 
Sey 4 I 


Taking exponentials of both sides and multiplying by y + 1 
gives 


y=(y+ loexp(x + C). 
Solving for y gives the general solution 


_ exp(x + C) 
—~{—exmp(e+C) 
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(exp (x + C))? 
(1 — exp (x + C)) 


(11) Using the partial fraction expansion found in 5S(ii) the 
equation can be written 


a ‘ 1 
dx 3(x-—1) 3(x+2) 
Integrating both sides gives the general solution 


y = 4log, (1 — x) + flog, (x +2)+C 
(since x + 2>0and 1 —x> 0). 


dy _ exp (x + C) 


Check: — = 
dx 1-—exp(x+C) 


Check: 
i pe as 
dx 3(x-—1) 3(x+2) 
x+1 
ee 


(iii) After dividing both sides of the equation by B* — y’, 
we can integrate both sides to obtain 


A 
jo zdy = ax + C. 


Using = partial fraction expansion obtained in Exercise 
5(iv) this can be written 


1 1 
Sag tama)?" Jace, 


Carrying out the integrations we obtain 


A 
—|] B = & 

~ Slog, (B—y)+ ap Oke (B+y)=x+ 

(since B—y>QOand B+y>Q0). 

This can be written 


By ie 


log. = ae + C) 

ne) (B + y) = (B— y)exp (2B(x + C)/A) 
Bexp (2B(x + C)/A)— B 

hence — 


1 + exp (2B(x + C)/A) © 


d 
Check: calculation of both and B? — y* gives 


B? exp (2B(x + C)/A) 
(1 + exp (2B(x + C)/A))* 


Solution to the exercise in Section 4 
1. (i) This equation can be written in the form 


dy d 
poy + PO) = ax) 


= 1. It 


d 
by choosing p(x) = x and q(x) = 1, for then oe 
follows that the equation can be written 


d(xy) 
dx 


and can therefore be integrated directly to give 


=i. 


xy=x+C. 


(ii) Here p(x) = sin x, Sad = cos x, and q(x) = sinx and 
a 
the given differential equation can be written 
d((sinx)y) 
——_——— =sinx. 
dx 


It can therefore be integrated directly to give 


(sinx)y = —cosx + C. 


=yty’. 
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(iii) Here p(x) = e* ecihcd, = 3e°* and q(x) = e** and the 
x 
given differential equation can be written 
dey) ox 
dx 


It can therefore be integrated directly to give 


er*y =te* +C. 


Solutions to the exercises in Section 5 


1. (i) A, (methods B and C could also be used but A is the 
most straightforward), (ii) B or C, (i) B or C, (iv) C, (v) B, 
(vi) D, (vit) C, (viii) B, (ix) D. 


2. 


3. (i) Euler’s method gives the recurrence relation 
41 =J),+01 x (1 —x,y,) 
where x, = r x 0.1. The method starts from x, = 0, yy = 1. 


(ii) the calculation is given in the following table 


1.1 
1.189 
1.26522 


11 : 89 ee 
0.2378} 0.7622 0.07622 
3 0. 3} 1.26522 


The Euler’s approximation to y(0.3) is shown in bold-face 
type. 

4. (i) Since the error is roughly proportional to the step 
length, the error would be approximately halved. 


(1) The number of calculations would be approximately 
doubled. 


6. Substituting the condition y(4) = 8 into the general 
solution gives 4C + 16 = 8 therefore C = —2 and so the 
required particular solution is 


y = —2x + x’. 

7. (i) The standard form of this equation is 
ay | 
ie 
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Using the separation of variables method gives 


1 
—— = — (since y > 0) 


thus [ja-- | dx + C. 
hence ——= —log,x+C (since x > 0). 
y 


The general solution is therefore 


1 
ae log,x —C 
dy —1 y? 
Check: — = —————___ = ——. 
en ax ~ x(log, x —C)? < 
(11) The standard form of this equation is 
EME 
dx x 


Since this is linear we can use the integrating factor method. 


The integrating factor is 
2 
exp [2 ix = enp(2iog x)= x*. 
x 
Multiplying our original equation by this factor gives 


so 2xy =x" 
dx 
which can be written 
axe =. 
= x”, 
dx 


Direct integration gives x7y = 4x* + C and so the general 
solution is 


“3 
y= 4x? +o. 
ers +2 - aes eres 
ea 2 2 


(111) This equation can be solved using separation of 


variables. Since we are only looking for positive solutions we 


can divide by y + y” and integrate to obtain 


[os : sdy = [dx +. 
a 3 


The right-hand side of this equation can be integrated to give 


x + C. The integration on the left-hand side can be carried 
out using partial fractions: let 


ae ee 
04> Ly Se 
then 


=N,(i+y)+ N.y 
=N,+(N,+N>)y. 
It follows that N, = 1 and N, = —1. Hence 


1 1 1 
d = {(-— d 
Fess F es ! 


= log. y ae log, (1 + y) 


ales)) 
“\i+y 


Equating the two sides of the equation gives 


(since y > 0) 


y 
log, |] = x + C 
“(7] 


; y 
1.e. — =exp(x + C). 
ee p ( ) 
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The general solution can therefore be written 
exp + C) 
1 —exp(x+C) 
dy _  exp(x + C) 
dx 1-—exp(x+C) 


(exp (x + C))? 
(1 — exp (x + C))? 


(iv) This equation is linear and can be solved using the 
integrating factor method. The integrating factor is 


exp [ fi-24 ae Fi 


Writing the equation in the form 


Check: =yt+y’. 


dy 
A eg, OPE 
dx ‘ 
_ and multiplying by the integrating factor gives 
d 
on 2x ly = 2go Py or —e” **x 
dx 
which can be written 
ae ee 
dx 
SO e- ys + [e-rxax +C. 


Integrating by parts we obtain 

e *y =ixe * — [3-2 ax +C 
=4xe * 4+4e°2 +C. 

The general solution is therefore 

y=3x+1+ Ce. 

Check: we have 


also 
2y —x=x+5+4+2Ce* —x=14+4 2Ce* 


and both sides agree. 


8. The integration can be carried out using partial fractions. 


Since x* + 3x + 2 = (x + 1)(x + 2) we can write 
x+4 N, N, 


Ox oat i ox 42 


therefore 
X+4=N,(x+2)+N,(x + 1) 
=(N,+N,)x + (2N,+N,) 


hence 
N, +N,=1 


The solution of this simultaneous equation is N, = 3. 
N, = —2 and so 


x+4 i 3 Z 
oe ge 


(Check: if x = —4 then both sides of this equation vanish.) 
We can now carry out the integration 


x+4 3 2 
(~~, ax = | — ——|dx 
x* +3x+2 ta ON ae oe ee 


= 3log, (—(x + 1)) — 2log,(—(x + 2))+ C, | 


where we have used the fact that x < —2 implies (x + 1) and 
(x + 2) are negative. 
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Solutions to home exercises on RECREL in Section 6 


1. (i) The complete table is: 


Information 


Formula for m(x, y) X+Y 
Method (Euler) none required 
Initial value of y 0) 
Initial value of x 0 
Number of steps 5 

Step size 0.2 


Type of print-out none required 


(ii) The response to Option 36 becomes 0.002; Option 4 | is 
replaced by Option 40, and the response to it is 50 because 
we only want every 50th calculated value (since 50 x 0.002 

= 0.1, the interval between successive printed values of x). 


2. (i) A suitable table is 


Information Response 


Order of recurrence relation 
Recurrence relation is 
Method (forward recurrence) 
Initial value uy 

Number of terms 

Type of print-out 


1 —(R + 1)*U(R) 
none required 
0.63212056 

25 
none required 


(ii) Part (i) can be left in the machine but the following new 
options must be inserted. 


Information Response 


none required 
any number 
between 0 
and 1 say 
35 


Method (backward recurrence) 
Value of uy 


Number of terms 


Note: Inserting Option 21 automatically cancels Option 20 
given in part (i), also, the response 35 to Option 35 
automatically cancels the previous response of 25 to Option 
35 given in part (1). 
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